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A novel model of incompressible magnetohydrodynamic turbulence in the presence of a strong 
external magnetic field is proposed for explanation of recent numerical results. According to the 
proposed model, in the presence of the strong external magnetic field, incompressible magnetohy- 
drodynamic turbulence becomes nonlocal in the sense that low frequency modes cause decorrelation 
of interacting high frequency modes from the inertial interval. It is shown that the obtained non- 
local spectrum of the inertial range of incompressible magnetohydrodynamic turbulence represents 
an anisotropic analogue of Kraichnan's nonlocal spectrum of hydrodynamic turbulence. Based on 
the analysis performed in the framework of the weak coupling approximation, which represents one 
of the equivalent formulations of the direct interaction approximation, it is shown that incompress- 
ible magnetohydrodynamic turbulence could be both local and nonlocal and therefore anisotropic 
analogues of both the Kolmogorov and Kraichnan spectra are realizable in incompressible magne- 
tohydrodynamic turbulence. 
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I. INTRODUCTION 

Although magnetohydrodynamic (MHD) turbulence 
has been extensively studied for the last 40 years U [j , 
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many physical aspects of the problem still remain unclear 
(for recent reviews see Refs. [ll IM Ell III ) ■ The first 
model of incompressible MHD turbulence was proposed 
by Iroshnikov 1] and Kraichnan Q. The Iroshnikov- 
Kraichnan (IK) model of MHD turbulence is based on 
the so-called Alfven effect [lj| - nonlinear interaction is 
possible only among Alfven waves propagating in oppo- 
site directions along the mean magnetic field. Therefore, 
IK model assumes that the energy cascade in MHD tur- 
bulence occurs as a result of collisions between oppositely 
propagating Alfven waves. Consider isotropic excitation 
of Alfven waves on some outer scale with a character- 
istic velocity vo "C Va, where Va is the Alfven velocity. 
IK model assumes that the energy transfer is local and 
isotropic in the wave number space. The characteristic 
time scale of the Alfven wave collision is t^;? ~ (VaA;) -1 , 
where k is the wave number. Using the governing equa- 
tions of incompressible MHD it can be shown that during 
one collision distortion of each wave packet Svi is of order 
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Because these perturbations are summed with random 
phases N ~ (vi/Svi) 2 ~ (Va/vi) 2 collisions are necessary 
to achieve the distortion of order unity. Therefore for the 
energy cascade time t^ s we have 



dimensional energy spectrum, we obtain 
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which represents IK spectrum of incompressible MHD 
turbulence. 

IK model of the MHD turbulence is isotropic. How- 
ever, presence of a mean magnetic field has a strong ef- 
fect on the turbulence properties, in contrast to a mean 
flow in hydrodynamic turbulence, which can be elimi- 
nated by a Galilean transformation. The anisotropy of 
MHD turbulence had been seen in various numerical sim- 
ulations 0, a a HI El- A theory of anisotropic MHD 
turbulence was proposed by Goldreich and Sridhar (GS) 
|5| . GS model implies that the dynamics of turbulence is 
dominated by the perpendicular cascade (i.e., by the cas- 
cade in the q-space, where q is the component of the wave 
number vector perpendicular to the mean magnetic field) , 
whereas the parallel size of turbulent 'eddies' (wave pack- 
ets) is determined by critical balance condition which im- 
plies that the characteristic time scale of wave collision 
(pV^) -1 , where p is parallel wave number, is equal to 
the characteristic time scale the energy cascade T cas . GS 
model also assumes that the energy cascade is local in the 
wave number space and therefore r^f ~ (viq)^ 1 . This 
yields (we suppose that turbulence is generated at per- 
pendicular scales l qo ~ 1 /go with a characteristic velocity 
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which represents the anisotropic analogue of the Kol- 
mogorov spectrum of isotropic hydrodynamic turbulence. 
The anisotropy of fluctuations is determined by the criti- 
Taking into account the relations e ~ vf /r cas , vf ~ fc-Efe, cal balance principle, which yields p ~ q 11 , with fj, = 2/3, 
where e is the energy cascade rate and is the one and therefore as the cascade proceeds to larger q, the 
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eddies become more elongated along the direction of the 
mean magnetic field. 

however, GS model does not totally agree with the 
results of the recent numerical simulations [9|, [ill, [IH, . 
It has been found [9( that in the presence of the strong 
external magnetic field incompressible MHD turbulence 
is strongly anisotropic with // « 1/2 and has the one 
dimensional spectrum similar to IK one £7. ~ q~ 3 ^ 2 . As 
it was shown in Ref. [11] the anisotropic spectrum of 
incompressible MHD turbulence strongly depends on the 
strength of the external magnetic field. It was found that 
the perpendicular scaling of fluctuations changed from 
GS form to the form found in Ref. Q as the external 
field was increased from 7 = B 2 /pv 2 <C 1, where Bq is 
the external magnetic field and p is the density, to the 
limit 7 ^> 1. 

For resolution of this inconsistency it has been noted 
fl4| that if in the presence of the strong external mag- 
netic field nonlinear interaction of the counter propagat- 
ing fluctuations is reduced by the factor vi/Va, so that 
the perpendicular cascade timescale becomes 
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then for the one dimensional spectrum one obtains 
Ef IK ~ {eV A ) 1/2 q~ 3/2 - This model of MHD turbulence 
represents an anisotropic analogue of IK model. The 
characteristic parallel scale of the turbulent eddies is de- 
termined by the condition similar to GS critical balance 
condition l/p ~ Vat ccls , which yields p ~ q 1 / 2 . 

Although anisotropic IK model properly represents 
scaling indices of incompressible MHD turbulence in the 
presence of the strong external magnetic field, it still con- 
tradicts with some other features of the turbulence ob- 
served in the numerical simulations. Indeed, according 
to Eq. ([5]), the ratio of the cascade timescale to the 
eddy turnover time r to ~ (fc^z) -1 is T cas /T to ~ vi/Va- 
For the parameters used in Ref. [§] the ratio Vi/Va ~ 
10~ 2 — 10 _3 ,whereas the ratio r cas /r to ~ 1 (see Sec. 
5.6.3 of Ref Q). Therefore, anisotropic IK model gives 
strongly understated estimate for the characteristic cas- 
cade rate. 

In this paper we present a novel model of incompress- 
ible MHD turbulence. The main assumption of the model 
is the nonlocal character of MHD turbulence in the pres- 
ence of a strong external magnetic field. The model im- 
plies that low frequency modes cause decorrelation of in- 
teracting high frequency modes from the inertial inter- 
val. Proposed model represents an anisotropic analogue 
of Kraichnan's nonlocal spectrum of hydrodynamic tur- 
bulence [23l | and properly reproduces main features of in- 
compressible MHD turbulence observed in the numerical 
simulations. 

Our study of incompressible MHD turbulence is per- 
formed in the framework of Kadomtsev's weak coupling 
approximation (WCA) [24j , which represents one of the 
equivalent formulations of Kraichnan's direct interaction 
approximation (DIA) [23j. Our analysis does not al- 



low us to prove that in the presence of a strong ex- 
ternal magnetic field incompressible MHD turbulence 
becomes nonlocal. However, performed analysis shows 
that for the systems with more then one degrees of free- 
dom the concept of adiabatic and resonant interactions 
does not hold and consequently, in contrast to isotropic 
hydrodynamic turbulence which is local, incompressible 
MHD turbulence could be both local and nonlocal. Ob- 
tained local solution reproduces GS model of MHD tur- 
bulence, whereas obtained nonlocal solution represents 
an anisotropic analogue of Kraichnan's nonlocal spec- 
trum of hydrodynamic turbulence. 

The paper is organized as follows. The heuristics of 
the proposed model is presented in Sec. HU Physics of 
the decorrelation mechanism is discussed in Sec. IIIII The 
WCA equations of incompressible MHD turbulence are 
derived in Sec. [IVJ Analysis of the WCA equations of 
isotropic hydrodynamic turbulence is given in Sec. [V] 
Local as well as nonlocal solutions of the WCA equa- 
tions of incompressible MHD turbulence are obtained in 
Sec. I VII Locality of incompressible MHD turbulence is 
discussed in Sec. IVII1 Conclusions are given in Sec. |IXJ 



II. HEURISTICS OF THE MODEL 

Kraichnan obtained nonlocal spectrum of hydrody- 
namic turbulence as a possible inertial range solution of 
the Eulerian DIA equations [23]. Kraichnan's nonlocal 
spectrum of isotropic hydrodynamic turbulence can be 
derived based on the following qualitative arguments [25| : 
assume that the autocorrelation time scale is determined 
by advection of a small-scale eddies (from the inertial 
range) of the size fc _1 through a distance of the order of 
its own size by a macroscopic flow of rms velocity which is 
determined by the characteristic velocity of energy con- 
taining eddies Vq, i.e., assume that T* r ~ l/(v Q k). Then 
the energy transfer rate reduces by the factor T^ c r /T to 
compared to the Kolmogorov theory (here r to ~ l/(vik) 
is the characteristic eddy turnover timescale) and there- 
fore the characteristic timescale of the turbulent cascade 
becomes 
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For the velocity fluctuations and the one dimensional 
energy spectrum one obtains vi ~ I 1 / 4 and Ejf ~ 
(ewo) 1 ^ 2 ^ -3 ^ 2 , respectively. This represents Kraichnan's 
nonlocal spectrum of isotropic hydrodynamic turbulence 
& 

As it was proposed by Richardson [26| and confirmed 
by various experiments (see, e.g., Refs. |27l. [28|] and ref- 
erences therein) the energy cascade in isotropic hydrody- 
namic turbulence is local and therefore the Kolmogorov 
spectrum is the only one that is realized in the inertial 
range of isotropic hydrodynamic turbulence. 

We propose that in contrast to hydrodynamic turbu- 
lence, incompressible MHD turbulence is non-local in the 
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presence of the strong (7 > 1) external magnetic field, 
i.e., an influence of low frequency modes leads to decor- 
relation of interacting high frequency modes from the 
inertial interval. Physics of the decorreelation mecha- 
nism is discussed in the next section. We argue that 
E q ~ <7~ 3 / 2 one dimensional spectrum observed in the nu- 
merical simulations represents the anisotropic analogue 
of the Kraichnan nonlocal spectrum. In particular, sim- 
ilar to GS model, presented model implies that incom- 
pressible energy transfer is strongly anisotropic and dom- 
inated by perpendicular cascade. The later cascade is 
assumed to be nonlocal, and similar to Eq. ([5]) we have 
for the characteristic cascade timescale 



1 ^0 



(7) 



therefore for the velocity fluctuations and the one dimen- 
sional energy spectrum our model implies 
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The relation between the parallel and perpendicular sizes 
of turbulent eddies is determined based on the following 
arguments the lifetime of a wave packet is of order 
r cas and therefore for the correlation lengthscale along 
the field line we have 
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Using Eq. this yields 
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Presented model properly reproduces main character- 
istics of incompressible MHD turbulence observed in the 
numerical simulations [jj [Tl| . 

Explicit dependence of the cascade rates given by Eqs. 
© and on the characteristic velocity of low frequency 
energy containing modes represents important feature 
of Krichnan's nonlocal spectrum as well as considered 
model and underlines the nonlocal character of the tur- 
bulence. 



III. PHYSICS OF THE DECORRELATION 
MECHANISM 

The governing equations of incompressible magnetohy- 
drodynamics in the Elsasser variables are 

9 t U = -(W- V)U- VP + i?V 2 U + P m V 2 W, (11) 



d t W = -(U • V)W - VP + i/ ro V 2 U + DV 2 W, (12) 



where U = v + b and W = v — b are the Elsasser vari- 
ables, v is the velocity field, b = B/^iirp is the magnetic 
field in velocity units, P is the total (hydrodynamic plus 
magnetic) pressure normalized by the density, dt = d/dt, 
v = [y + v m )/2, T>ra = (v — u m )/2, v is the kinematic vis- 
cosity and v m is the magnetic diffusivity. 

Consider the packet of Alfven waves propagating along 
the uniform magnetic field bo || 2, with the characteris- 
tic magnetic field perturbation bo _L z, the velocity per- 
turbation vo = —bo and the characteristic parallel and 
perpendicular wave numbers po and qg, respectively. On 
this background consider a high frequency wave packet 
propagating in the opposite direction with the character- 
istic parallel and perpendicular wave numbers p u po 
and q u 3> qo, respectively. Denote the Elsasser variable 
associated with the high frequency wave packet by u„. 
Assuming bo as approximately constant and dropping 
dissipation terms Eq. (jlip yields 



(dt - b • V - b • V + v ■ V)u„ « 0, 



(14) 



Noting that vo = —bo, this equation gives for character- 
istic frequencies of the high frequency packet 



-b pu + 2v q u . 



(15) 



If we consider dynamics of a high frequency Alfven wave 
packet propagating parallel to z-axis (with the charac- 
teristic parallel and perpendicular wave numbers p w ~ 
p u ~ p and q w ~ q u <~ q respectively), then similar con- 
sideration yields for the frequency 



boPv 



(16) 



v-w = v-u = o, 



(13) 



Therefore, presence of the low frequency wave packet do 
not influence propagation of the high frequency packet 
along 2-axis (Alfven effect) whereas the packet propagat- 
ing in the opposite direction is moved with the velocity 
2vo in the direction perpendicular to z-axis. This cir- 
cumstance represents the bases for understanding, why 
incompressible MHD turbulence becomes nonlocal in the 
presence of the strong external magnetic field. Indeed, 
consider two interacting packets of high frequency alfven 
waves moving in opposite directions on the same back- 
ground as above (i.e., in the presence of strong exter- 
nal magnetic field and the low frequency Alfven wave 
packet) . The perpendicular length scale of the packets is 
of order 1 j q and consequently, due to the described above 
action of the low frequency packet, the interacting pack- 
ets would be pulled apart during the time scale l/(qv ). 
Therefore, the characteristic timescale of the unit act of 
interaction shortens compared to the Kolmogorov auto- 
correlation timescale l/(qvi). As it was shown in the 
previous section this automatically implies that the tur- 
bulence would have the spectral characteristics given by 
Eqs. 

In the consideration above we implicitly assumed that 
the low frequency packet do not contribute to the mean 
magnetic field which acts on the interacting high fre- 
quency wave packets (i.e., it has been assumed that the 
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anisotropic high frequency wave packets are formed along 
the external magnetic field b ). If this is not the case then 
the situation changes entirely. Indeed, assume that the 
low frequency wave packet contributes to the mean field 
acting on high frequency packets, i.e., assume that the 
high frequency packets are formed along the axis bo + bo 
instead of bo. In this case, instead of Eqs. (| 1 5 [1 - (fl6]) 
similar analysis yield 

LOu ~ —(bo + b a )p u + v ■ q„, (17) 

« {bo + b )p (18) 

The second terms in Eqs. (fT7 |) - (fT5|) describes the Dopplcr 
shift caused by the velocity field of the low frequency 
wave packet and could be removed by corresponding 
Galilean transformation. Consequently, in this case the 
only effect caused by the presence of the low frequency 
modes is the change of the mean field, and they do not 
cause the decorrelation of high frequency modes. This 
arguments lead us to the conclusion that only the low fre- 
quency modes which do not contribute to the mean field 
cause decorrelation of interacting high frequency packets. 

Consequently, presented model implies that in the 
presence of the strong external magnetic field bo 3> bo 
the low frequency energy containing modes do not con- 
tribute to the mean field which acts on the high frequency 
fluctuations in the inertial interval. As a result, decor- 
relation mechanism described above causes formation of 
the anisotropic analogue of the Kraichnan nonlocal spec- 
trum. On the other hand, in the absence of the strong 
external magnetic field, the mean field is formed by the 
low frequency (energy containing) modes, the nonlocal 
decorrelation mechanism is not at work, turbulence is 
local and therefore is described by GS model. 

IV. THE WCA EQUATIONS FOR 
INCOMPRESSIBLE MHD TURBULENCE 

Consider incompressible MHD turbulence in the pres- 
ence of the constant magnetic field B directed along 
z axis. Presenting the Elsasser variables as a sum of 
the mean and fluctuating parts W = W + Wj, U = 
U + Ui, where U = and W = — V^, and 
= B /\/47rp is the Alfven velocity, and dropping 
dissipation terms Eqs. (|TT|) - (|T5]) yield 

d t Vx - V A d z V l = -(Wj • V)Ux - Vp (19) 

dtWi + Va<9 z Wi = -(Ui • V)Wi - Vp (20) 
Performing the Fourier transform defined as 

Uk uj = * J exp(iujt - ik ■ x)Ui (x, t)d 3 xdt : (21) 

and eliminating pressure terms we obtain 

(^ + ^ k )u k = Aui-fc(fc-ui)](k.w a )djf )2) (22) 



(^-^ k )w k = /[wi-£(£.wi)](k.u a )dFf i2 , (23) 

where k = (k, cj), the caret denotes the unit vector, Ui 
denotes u ki , Wk = Va^-z is the frequency of the Alfven 
wave, dT\z = d 4 k\d 4 k 2 <5 k _ kl _ k2 , and <5 k _ kl _ k2 = 
ki — k 2 ) is the Dirac delta function. 

Incompressible MHD turbulence is governed by inter- 
action of shear Alfven waves, whereas pseudo Alfven 
waves play a passive role [!, Q- Therefore, in the pre- 
sented paper we consider the shear Alfvenic turbulence. 
Defining the unit polarization vector of the shear Alfven 
waves e k = k x z, and introducing the amplitudes of the 
shear Alfven waves as 

w k = %e k , u k = ^k 6 k, (24) 
Eqs. (l22 |) -([23 |) reduce to the following equations 

/oo 
11,201^2^2, (25) 
-oo 

/oo 
21,2^1^2^,2, (26) 
-oo 

where Ti,2 = i(e k • e kl )(k • e k2 ) is the matrix element of 
interaction. 

To achieve any progress in analysis of Eqs. J25p - (l2l)|) 
some closure scheme should be used. In Ref. [5J the eddy 
damped quasi normal Markovian (EDQNM) approxima- 
tion was used for the analysis of Eqs. (f25|) - ([26|) . In the 
framework of this approximation the linear damping term 
is added to the equation for the third order moments of 
the turbulent fields and afterwards some assumptions are 
made regarding the eddy damping rate. But an assump- 
tion of some a priori form of the eddy damping rate auto- 
matically fixes the property of the locality of turbulence - 
is it assumed to be local or not. Note that locality of the 
turbulence is one of the main assumptions of GS model 
0]. Here we use the WCA (DIA) for the study of locality 
of the incompressible MHD turbulence. 

In the framework of the DIA [23| one differentiates di- 
rect nonlinear interactions among three Fourier modes 
with the complanar wave vectors (k + ki + k2 = 0) and 
indirect interactions when the modes are interacting by 
means of other Fourier modes. The idea of the closure 
of equations for the turbulent fields is based on the as- 
sumption that the influence of the indirect interactions 
on the turbulence dynamics can be neglected in compar- 
ison with the direct interactions. 

We use the standard WCA technique [24| for the 
derivation of the WCA equations of incompressible MHD 
turbulence. One of the main nonlinear effects described 
by the nonlinear interaction terms on the right hand sides 
of Eqs. (|2"5|) - (f2"6"|) is the nonlinear decay of the mode </> k . 
The intensity of this process is proportional to the am- 
plitude of the mode. The WCA implies the isolation 
of this part of the nonlinear interaction. For this pur- 
poses one should add to both sides of Eqs. (|2"5"|) - (f2"6")l the 
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terms iC^k an d "i-C^^Pki respectively. Multiplying ob- 
tained equations respectively by and ip^, (here and 
hereafter asterisk denotes the complex conjugated value), 
ensemble averaging, considering different wave modes as 
statistically independent and introducing the following 
notations 



we obtain 



J^ k _ k „ (27) 



ri, 2 <0i^ 2 C>^i fc 2, (28) 



T^iM&dF^. (29) 



The second terms on the right hand sides of Eqs. l|28jl - 
contain contributions from both direct and indirect 
nonlinear interactions. Following the procedure devel- 
oped in Ref. [24| for the elimination of the contribution 
of the indirect nonlinear interactions, let us represent the 
turbulent fields as 



k ' 



4" 



tS.o) 



(30) 



(31) 



where in the zeroth approximation the fluctuations are 
assumed to be uncorrelated. Then the third order mo- 
ments on right hand sides of Eqs. P8|) - (f!?9")) vanish in the 
zeroth order approximation. In the first order approxi- 
mation we have 



For 



Eq. ||25j gives 



( W -u, k )4 1)= / TwWWdfif. (33) 

To isolate the contribution of only direct nonlinear in- 
teractions we note, that the duration of the unit act of 
the nonlinear interaction is ~ VCg- Therefore, to 
hold contribution of only direct nonlinear interactions we 
should replace (ui — w k ) by (w - Wk + *C^) hi Eq. (|33|) . 
This yields 



(32) 



(°)„/,(°), 



1 



Similarly for we obtain 



/oo 
r ll2 0i o) 4 o) d^ ia . (34) 
-oo 



(1) 



1 



(w + w k + iC- ) J- 



T^^dF* (35) 



Substituting Eqs. (p ]) . ([34 ]) and ([35" ]) into Eqs. (l28l) - 
(|29| . considering the fields efc and i/'k as uncorrelated 
(i.e., assuming ((f>iip2) ~ 5], which physically implies 
that we consider MHD turbulence with zero residual en- 
ergy , or equivalently with equal kinetic and mag- 
netic energies), using for the forth order moments the 
Gaussian relations 



S^M^) = (<f>3<f>k) ("04-02) = IqI 2 * 



k-k 3 k 2 +k 4 



(36) 



A ^2 ^ki-k 3 ^k 2 -k 4 ' 



and taking into account that according to the definition 
of Ct, the right hand sides of Eqs. (|25 ]) -(p9" |) should not 
contain the terms proportional to ig we finally arrive at 
the following equations 



T1.2Tk.-2 



-00 u x - w kl + iCi" 

Jl,2^k,-2 



4" <tfT, 2 > 



- 1 2 u " r 1,2) 



(37) 



(38) 



(w - uik + it 



(UJ + LUk + i& 



+ \|2 



+ uikt + i(i 

/OO 
|T 1)2 | 2 7+/ 2 -^ )2 , (39) 
-OO 

/OO 
|T 1)2 | 2 7r/+rf^ 2 . (40) 
-OO 

Eqs. I[3"?| ) -(fl0"] ) represent the WCA (equivalently the 
Eulcrian DIA) equations for incompressible MHD turbu- 
lence. Noting that T\ 2 = T k _ 2 and defining 



w T ^k + iC 



± ■ 



(41) 



after straightforward manipulations Eqs. (|37p - ([40|) can 
be rewritten as 



i(w - Wk)rj±; 



/OO 
|Ti, 2 | 2 r+/ 2 -AFfc 2) (42) 
-OO 

/CO 
|T lj2 | 2 r r /+^ 2 , (43) 
-00 

/CO 
|T li2 | a J+J^dJ* )a - 
-OO 

/CO 
|T li2 | 2 r+I 2 -(^ 2 , (44) 
-00 

/oo 
-OO 

/OO 
|r lia | 2 rrtfd?? i2 , (45) 
-00 

The aim of the further analysis is the study of Eqs. 
([37 | - p0|) . or equivalently, (|42|) - (|45"|) and similar equations 
for isotropic hydrodynamic turbulence. But before per- 
forming this analysis we shortly discuss several topics 
related to the nature of the obtained equations. 
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A. The equivalence of the WCA and the DIA 

As it was mentioned above, the WCA is one of the 
equivalent forms of the DIA. Essentially the WCA is the 
DIA formulated in the frequency domain instead of the 
time domain (2f| [29[. To perform the DIA analysis of 
incompressible MHD turbulence one should start from 
Eqs. (J25j)-([26l) in the time domain (Eqs. (15) of Ref. 
Q ) , add random forcing terms and define the Green func- 
tions G + (k, t — t') and G~(k, t — t') of the obtained equa- 
tions. Then using standard assumptions of this approach 
(see, e.g., [25[ and references therein) for separation of 
the direct and indirect interactions one should obtain 
the closed set of equations for functions Q + {k, t,t') and 
Q~(k, t,t') defined as 

(<Mi)M*')> =Q + (k,M')<5k-k', (46) 



these equations allows us to develop the method for anal- 
ysis of the WCA equations for the simplest example - the 
inertial range of isotropic hydrodynamic turbulence. Sec- 
ondly, we derive Kraichnan's nonlocal spectrum. Com- 
parison of the results obtained in this section with the 
results for inertial range of incompressible MHD turbu- 
lence obtained in the next sections allow us to show that 
the energy spectrum E q ~ q~ 3 / 2 observed in the numer- 
ical simulations of incompressible MHD turbulence in a 
strong external magnetic field represents the anisotropic 
analogue of Kraichnan's nonlocal spectrum. 

The equations analogous to ([57]) and (|4"4")l for the iner- 
tial range of the hydrodynamic turbulence have the form 

J 5", rfe W,2- (48) 
Wi + wk{ + iQ 



(47) 



and the Green functions G ± (k, t-t'). In Eq. (gBJ) <f> k (t) 
is the inverse Fourier transform of 0^ with respect to 
u>. For stationary turbulence Q ± (k, t, t') — Q ± (k, t — t'). 
Close analogy between the WCA and the DIA could be 
seen if one notes that according to Eqs. ([27]) and (|4"6")h 1^ 



are the Fourier transforms of Q ± (k, t 
to t — t', where ' 
G ± (k,t-t') [1 



t') with respect 



whereas r?/27r are the Fourier transforms of 



B. The conservation laws 

It can be checked that Eqs. (|4"2")) - (r4"5"l) conserves total 
energy of the both types of the shear Alfven waves. In- 
deed, performing the inverse Fourier transform of Eqs. 
(|4"4")) - (|45p with respect to u>, setting the temporal vari- 
able equal to zero, integrating both sides of the obtained 
equations over the whole k space and taking into ac- 
count that 1^ = it can be shown that considered 
nonlinear interactions conserve total energy of the both 
types of the shear Alfven waves H + = 1/2 J |</)| 2 d 4 &; and 
H~ = 1/2 J |V>| 2 d 4 fc, or equivalently, Eqs. (|3I |) -(|l5" |) con- 
serve both the total energy H = H + + H~ and the cross 
helicity H = H+ - H~ . 



V. THE WCA FOR ISOTROPIC 
HYDRODYNAMIC TURBULENCE 



i {lo + ivk 2 ) 4 = fc 2 rjf 



/oo 
b^T^df^, (49) 
-oo 



where 



«1,2 = 2 



o (k.k0 2 (k.k 2 ) 2 

1.2 p 1.2 + 

(k-ki)(k-k 2 )(ki -k 2 ) 



h 2 h 2 h 2 



(k-k!) 3 (k-k 2 )(ki-k 2 ) 



k 4 k 2 



k 2 k 2 



(50) 



(51) 



and are related by the relation 141 [1 with uj k = 
—ivk 2 , and is denned by (vg ■ v^,) = I^S k _ k ,, where 
is the Fourier transform of the turbulent velocity field. 
Eqs. ([48 |) - ([5T|) are useless until some assumptions are 
made about the frequency dependence of and l£. 
Equivalently, in the framework of the DIA one should 
make some assumptions about the time dependence of 
G h (k, t - t') and Q h (k,t - t') (27j], which are the corre- 
sponding inverse Fourier transforms with respect to lo. 
One of the simplest and frequently used assumptions im- 

ply 0, S3 

G h (k, t-t') = exp [-|t£|(t - t')]H{t - t'), (52) 



Before performing analysis of the WCA equations (142]) - 
(|4"5)l for incompressible MHD turbulence, in this section 
we study much more simple problem - the WCA equa- 
tions for isotropic hydrodynamic turbulence. There are 
two reasons to perform this study. Firstly, although 
WCA equations for isotropic hydrodynamic turbulence 
have been derived by different authors 0, [3(3] , no meth- 
ods of the analysis have been indicated. The study of 



Q h (k, t-t')= exp [-|#|(t - t')]£ k , (53) 

where H(t) is the Heaviside (step) function, and £ k is 
the energy spectrum. Usually it is also assumed [U Hl| 
that the temporal autocorrelation scales of the spectral 
function and the Green function are equal 1 /r)£ = l/£ k = 
t^ c . From Eq. ([53")) it follows that temporal derivative 
of Q h (k, t — t') has the discontinuity at r = t — t' = 0. 
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As it was shown in Ref. [3l| this makes no problems in 
analysis of the turbulence dynamics if one assumes 



lim -£-Q h (k,r) =0. 
t^o dr 



(54) 



Similar to Eqs. ([52"]) - ([53"]) . in the case under considera- 
tion we assume 



7T UJ 2 + (rfc) 



h\2 ' 



(55) 



(56) 



Substituting these equations into Eq. I[i8"]) and integrat- 
ing over the frequency variables we obtain 



4 



ik 2 b lt2 £k 2 



d/Cf 2 , 



(57) 



where d/Cf 2 = d 3 kid 3 k 2 <$k-ki-k 2 - As we see the right 
hand side of Eq. (|57|) depends on w, whereas accord- 
ing to our assumption the left hand side does not. This 
means, that Eqs. ([55]) and ([56]) can not be valid for the 
whole (k, uj) space. The solution of this inconsistency is 
based on the following physical arguments: introducing 
the autocorrelation time scale r ac = 1 /rfc one also should 
assume that r ac is the shortest time scale in the problem 
and the dynamics of turbulence is totally determined by 
r ac . Alternatively, it can be shown that the main contri- 
bution in the energy balance equation (j4"9"]) comes from 
the modes with uj < r^, and therefore when studying the 
scaling properties of the inertial range of the turbulence 
one can set uj = in Eq. ([57]) . Taking also into account 
the isotropy of the turbulence (ry^ = 7/£ and ££ = we 
obtain 



4 



k~ 



-d/C* 



(58) 



It should be noted, that the standard DIA technique 
(see, e.g., Ref. [13|) also leads to Eq. ([58]) without ex- 
plicit assumption that r ac is the only timescale which 
totally determines the dynamics of turbulence in the iner- 
tial interval. As it is shown in the Appendix, the method 
of derivation of Eq. ([58]) in the framework of the DIA 
also implies implicitly that the formulated assumption is 
hold. 

For analysis of Eq. (|49p we perform inverse Fourier 
transform with respect to uj. In the right hand side of 
the obtained equation we can set v — 0, due to the fact 
that the dissipation range has negligible influence on the 
nonlinear energy transfer in the inertial range. Substitut- 
ing Eqs. ([55]) - ([56"]) . performing integration with respect 
to frequencies, using Eq. ([54"]) and taking into account 
the identity [13] 



1 



«i,2 = ~ (*>i,2 + h,i) ■ 



(59) 



we finally obtain 

vk 2 £ k = k 2 I b l:2 



>C_ 1.2 I . ^fc 2 (gfc! gfc) nyfc 



(60) 



For further simplification of this equation we multiply it 
by Airk 2 and integrate from some k in the inertial inter- 
val to the infinity. For isotropic hydrodynamic turbu- 
lence the effective dissipation takes place for high wave 
numbers and therefore 

/•OO /"OO 

Airv / k 4 £ k dk^u k 2 £ k d 3 k = £ h , (61) 

J k J —oo 



where Eh is the energy dissipation rate. Then Eq. I]60[) 
reduces to 

e h = A, I™ k^dk r h ffl\~ £ h k) dlCt 2 . (62) 

Jk J-oo »7fc+<+< 

To study the scaling properties of the turbulence we 
seek for the inertial range solution of Eqs. ([55]) and ([6"0]) 
in the form 



4 



B h k n . 



(63) 



(64) 



Substituting these expressions into Eq. I|58|) and taking 
into account that for k 3> k 2 , &1.2 ~ sin 2 a 2 , where a 2 is 
the angle between k and k 2 , we see that the integral on 
the right hand side of Eq. ([55]) diverges at k 2 — > if m < 
—3 and n > 0. Note that for the Kolmogorov spectrum 
to = —11/3 and n = 2/3. Therefore, the integral in 
Eq. (|58p is dominated by the low wave number energy 
containing modes for which Eq. (|63[) is not valid. 

To calculate the contribution of the energy containing 
modes we substitute Eq. ([6"3"]) into Eq. ([55]) and use delta 
function to integrate over k\. Noting that |k — k 2 | ~ k 
this yields 

B h k n » Bl x k 2 - n J £ k2 sin 2 a 2 d 3 k 2 ~ B^k 2 - n v 2 , 

where vq is the characteristic velocity of energy contain- 
ing vortices. From Eq. (]65p we obtain 



n = 1, 



B h ~ Vq 



(66) 



In contrast to Eq. {oSJ, due to the presence of the 
multiplier (£ kl —£ k ), integrals in Eq. ([62]) are convergent 
for the small values of k 2 when to > —4. Consequently, 
noting that the main contribution in the integral comes 
from the wave numbers with k\ ~ k 2 ~ k we obtain 



(67) 



and therefore 8 + 2to — n = and Ah ~ (ehBh) 1 ^ 2 . Using 
Eqs. (|66p and introducing the one dimensional energy 
spectrum E k = 4irk 2 £ k we obtain 



E k ~ ( £/l « ) 1/2 fc- 3/2 , Vk ~ kv , 



(68) 
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which represents Kraichnan's nonlocal spectrum [23| |. 
Note that in the considered model the spectral energy 
transfer is local [i.e., integrals in Eq. (|62|) are conver- 
gent] and the nonlocal character of the turbulence is re- 
lated to effective decorrelation caused by low frequency 
modes from energy containing interval [divergence of the 
integrals in Eq. (|58p and as a consequence shortening of 
the autocorrelation time scale from l/(kvi) to l/(kvo)]- 

Solution of the contradiction between Eq. (|68|) and 
Kolmogorov spectrum led Kraichnan to the formulation 
of so-called Lagrangian DIA 32j ■ Another approach was 
developed by Kadomtsev [24 1. He argued that analysis 
of Eq. ([58)) which led to Eq. (|65|) and ([66]) overestimates 
the contribution of nonlocal interactions. It was asserted 
that when studying nonlinear interactions of the modes, 
one should distinguish two kinds of nonlinear interactions 
- resonant and adiabatic. In Eq. (|58[) it is implicitly 
assumed that all the nonlinear interactions are resonant. 
In fact, interactions of the modes with very different wave 
numbers are adiabatic rather then resonant: the influence 
of the large scale fluctuations on the low frequency ones 
leads to adiabatic (WKB) change of the frequency u> and 
the wave number k of the high frequency fluctuation, 
similar to the weak inhomogeneity of the mean fields. 

As it was shown in Ref . [2J] , the dynamics of the tur- 
bulence in the inertial range is mainly determined by the 
resonant nonlinear interactions whereas influence of the 
adiabatic interactions can be neglected. From mathemat- 
ical point of view this means, that instead of integration 
over the whole k-space in Eq. one should restrict 

the area of integration by the area ak < | ki,2 1 < k/a, 
where a ~ 1, say a = 1/3, H or a = l/s/2, |]J,[I3. If 
so, using Eqs. ([63l) - ([64|l . instead of Eq. ([65]) we obtain 



Bk n ~ AB- l k b+r ' 



(69) 



2 /3 /fc- 5 / 3 . Note that for the 



Combining this equation with Eq. (|67[) we arrive at the 
Kolmogorov spectrum £7. 
Kolmogorov spectrum t^ c ~ T^ as ~ 1/kvi. 

As it will be shown in the Sec. IVII1 although the con- 
cept of resonant and adiabatic interactions yields right 
result for inertial range of isotropic hydrodynamic tur- 
bulence, for the systems with more then one degrees of 
freedom it does not hold. 



VI. NONLOCAL SPECTRUM OF THE 
INCOMPRESSIBLE MHD TURBULENCE 



7T ( WTWk ) 2 + (£k) 2 ' 



(73) 



In the presented paper we consider the symmetric case 

vi=^ = Vu, £^ = £^=£ k , (74) 

which corresponds to the turbulence with zero cross he- 
licity (££ — £ k = 0). In this case l£ = I^_ u , and 
it can be readily checked that Eqs. ([3"7| and ([3"9")l coin- 
cides with Eqs. and ([4"0)) respectively. Therefore, in 
the symmetric case there remain only two independent 
equations. 

Here we consider strong turbulence, i.e., assume that 
the autocorrelation timescale r ac = is the shortest 
timescale presented in the problem which totally deter- 
mines the dynamics of the turbulence. Taking this into 
account and performing integrations with respect to fre- 
quencies in Eqs. (|37p and (|44[) similar to hydrodynamic 
turbulence we obtain 



iri. 



vk 2 £\. 



\rp |2 ^k 2 (^ki £k) Alr k 
\Tx,2\ ; : d/C 12 - 



(75) 



(76) 



For analysis of anisotropic turbulence we assume that the 
dynamics is dominated by cascade perpendicular with re- 
spect to the mean magnetic field and therefore we assume 
r?k = rj q . For energy spectrum we assume Q 



~ A 



-"/( — 

1 M» 



(77) 



where f(u) is a positive symmetric function of u which 
becomes negligibly small for u 1. For u < 1, f(u) ~ 
1 such that f^ ^ f(u)du = 1. We also define the two 
dimensional energy density as 



£ k d P = At 



(78) 



Performing the integration with respect to pi and P2 in 
Eqs. ([75] ) - ([76] ) . multiplying both sides of Eq. ([76]) by 2?rg 
and integrating over p and q from some value in the iner- 
tial range to infinity and defining the energy dissipation 
rate 



For further analysis of Eqs. ([3T)) - (r4T)]) . or equivalently 
(|4"2")) - (|45p . similar to the hydrodynamic case we assume 

G ± (k,r) =exp(-|r ? ±|r±ia; k T)fl'(T), (70) 

g ± (k,r) =cxp(-\^\r±iuj k T)£ k , (71) 
or equivalently 

d=vt (72) 



we obtain 



>h 



e = 2ir I qdq 
'1 



fc 2 f k d 3 k = 2-kD 



\Tx,: 



Vqi + Vq: 



1 2 £<12 (£qi 



q 3 £~ddq. 



£a 



(79) 



(80) 



Vq + Vqi + Vq 2 



dQf, 2 , (81) 
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where dQ\ 2 = d 2 qid 2 q 2 <5 q _ qi _ q2 . 

Similar to Eq. we assume r\ q = Bq n . Taking 

into account that | !Z\ , 2 1 2 = Q 2 cos 2 0\ sin 2 6* 2 , where 0i i2 
are the angles between q and qi. 2 respectively, we obtain 
that the integral on the right hand side of Eq. (|80f is 
divergent for g 2 — ► if m < —2. Further analysis is 
analogous to one performed in Sec. [V]for hydrodynamic 
turbulence. Similarly we obtain: 



1/2-5/2 



(82) 



For the one dimensional spectrum E q = 2nq£q this 
yields E q ~ (evoY^q^ 3 / 2 . For the characteristic cascade 
timescale we obtain: 



1 vq 
qv q v q ' 



(83) 



The analysis of the anisotropic turbulence is not full 
until some relation between parallel and perpendicular 
length scales of the turbulent wave packets are specified. 
Technically this implies determination of \i and A. This 
could be done based on arguments discussed in Sec. [Til 
and lead us to Eq. (JTU|) or equivalently 



/! = 1/2, 



a "0 1/2 

A =^o ■ 



(84) 



VII. RANDOM GALILEAN INVARIANCE AND 
LOCALITY OF INCOMPRESSIBLE MHD 
TURBULENCE 

The study of the locality of incompressible MHD tur- 
bulence performed in this section is based on the methods 
developed in Rcf. 24]. Let us isolate resonant inter- 
actions in Eq. (|25|) and represent this equation in the 
following form 

(w-a; k )^k= / T\2$-l<^2&F\ ,2 + / ^l^i^d-F^ + 
Jn Jv 

[ Ti jk _i0iVk-id s ki + / T k _ 2 , 2 (/. k „ 2 ^ 2 d 3 kd;85) 

J Ax JA 2 

In this equation 1Z denotes the resonant area where 
qi ~ qi ~ q (as in the previous section we also assume 
that there exist some relation between the perpendicular 
and parallel length scales of turbulent wave packets [f|), 
T> denotes the high frequency area where q\ ~ g 2 ~> q, 
and A1.2 denote the areas qi.2 <C q, respectively. The 
analysis performed in the previous section shows that 
the area D has negligible contribution to the integrals 
in Eqs. ([75]) - ([76]) and will be neglected in the further 
consideration. However, it should be noted, that nonlo- 
cal interactions with small scale fluctuations could have 
important contributions in the dynamics of MHD turbu- 
lence with nonzero helicity p^Ult^. 

Performing Taylor expansion in the integrands of the 
third and forth terms on the right hand side of Eq. (1551) 
with respect to the small parameter (71,2/9 respectively 



and taking into account the definition of Ti. 2 and Eq. 
we obtain 

(lu - Wfc)^ = k • (Vj, + B L )4>y. + G00 + Gfrp + 

Ti, 2 0lV2<tf?2- ( 86 ) 

Similar manipulations of Eq. (|26[) yields 

(uj + Wk)ip£ = k • (V L - B L )i/> k + + G^ + 



Jn 



(87) 



In Eqs. ([5rJ| and ([57]) Gij are the First order terms of 
the Taylor expansion G^ ~ deft/dk and G^ ~ dip/dk. 
V l = J kl<k v kl ci 4 k x and B L r fci<fc bk lC i 4 ki are the ve- 
locity and the magnetic fields of the low frequency modes. 

According to the principle of resonant and adiabatic in- 
teractions, first terms of the right hand sides of Eqs. ([55)1 
and (|87p describe the transfer of high frequency pertur- 
bations by the low frequency modes and one should elim- 
inate these terms before applying the closure procedure 
described in Sec. [TV] In the case under consideration 
this could be done as follows: the terms proportional to 
Vj, can be eliminated by means of the Galilean transfor- 
mation (this will lead to the frequency renormalization 
lu' = uj — k ■ Vi), whereas the terms proportional to 
Bi can be included to w k . If so, the influence of low 
frequency modes on the dynamics of the high frequency 
modes comes to slow (adiabatic) deformation of the high 
frequency modes, described by G^ terms in Eqs. ([55]) 
and (|87[) . In the zeroth order approximation the adi- 
abatic interactions can be neglected and therefore Eqs. 
(|55 ]) -([87 1) formally coincide with Eqs. J55J-J2B]), except 
the integration now is performed only over the resonant 
area 1Z. Applying the WCA (DIA) closure scheme de- 
scribed in Sec. [TV] to these equations, we obtain Eqs. 
(|4"2")) - (|4"51) with the integration performed only over the 
resonant area 1Z. Then analysis similar to the performed 
in the previous section would lead to the GS model of 
the turbulence described by Eqs. Q. 

The WCA (Eulerian DIA) suggests that low frequency 
fluctuations always cause decorrelation of interactions of 
high frequency fluctuations and as it is known this result 
is wrong [24], |32| . In contrary, the principle of resonant 
and adiabatic interactions states that the influence of 
low frequency modes on the dynamics of high frequency 
modes always have adiabatic character. Although this 
principle yields right results for inertial range of isotropic 
hydrodynamic turbulence as well as for any system with 
one degree of freedom [i.e., when there exist one starting 
equation similar to Eqs. (|25 p - ([26]) ]. in general case this 
principle is not correct. 

Indeed, in the case of the velocity field the necessity of 
elimination of the terms proportional to V^ from Eqs. 
([86]) - ([87]l is caused by the requirement of the invariancc 
of the governing equations with respect to the random 
Galilean transformation [32] . According to this principle 
the 'right' form of Eqs. (|5BJ)-|(57J) has V L = 0. In the case 
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of magnetic field there exists no such a fundamental prin- 
ciple and consequently the necessity of the elimination of 
the term proportional to Bl is not obvious. Physical 
picture is the same as it was discussed in Sec. IIIII If the 
low frequency energy containing modes contribute to the 
mean magnetic field acting on the high frequency modes, 
then the terms proportional to should be included in 
terms proportional to lo^ on the left hand side of Eqs. 
(I8"6"|) - (j87l) . In this case the DIA closure scheme presented 
in Sec. IIVI and analysis similar to one presented in Sec. 
I VII leads to the GS model of the incompressible MHD tur- 
bulence. On the other hand, if the low frequency modes 
do not contribute to the mean magnetic field, then the 
same kind of analysis leads to the anisotropic analogue of 
Kraichnan's nonlocal spectrum. Therefore, we conclude 
that in contrary to the incompressible hydrodynamic tur- 
bulence, incompressible MHD turbulence could be both 
local and nonlocal. 

In the general case it is clear that the principle of reso- 
nant and adiabatic interactions holds for any system with 
one degree of freedom, since in this case the zeroth order 
term of the Taylor expansion describing influence of the 
low frequency modes on the high frequency modes always 
can be eliminated by corresponding Galilean transforma- 
tion. In the case of the system with two or more degrees 
of freedom random Galilean invariance can not guarantee 
the elimination of all zeroth order terms, as in the con- 
sidered above case of Eqs. (|86|) - (|87| and therefore the 
low frequency modes can cause effective dccorrclation of 
high frequency fluctuations. 



VIII. TWO DIMENSIONAL MHD 
TURBULENCE 

It is well known that dynamics of two dimensional hy- 
drodynamic turbulence is basically different from three 
dimensional turbulence. The origin of this difference 
is that in two dimensional case in addition to the en- 
ergy there exist another conserved quantity - enstro- 
phy. The situation is different in magnetohydrodynam- 
ics. As it was discussed above, MHD turbulence is highly 
anisotropic, dominated by turbulent cascade in the direc- 
tion perpendicular to the mean magnetic field. Therefore 
it has been argued that the local two dimensional dynam- 
ics in the plane perpendicular to the mean magnetic field 
describe an essential part of the three dimensional tur- 
bulence (especially in the presence of the strong mean 
field) and could be studied by two dimensional MHD. In 
this context two dimensional MHD turbulence has been 
extensively studies by different authors [33j, [3J, l35|. It 
seems interesting to compare presented model with the 
features of two dimensional MHD turbulence observed in 
numerical simulations. Numerical simulations performed 
in Refs. 33, 35] show that inertial range spectrum of two 
dimensional MHD turbulence in the absence of a mean 



magnetic field has the form 

E k = C^Etfk-*'*, (88) 

with Cfc ~ 1.8 and Em is the density of the magnetic 
energy. Noting that Em ~ Vq we conclude that this result 
coincides with the prediction of the presented model for 
perpendicular cascade given by Eq. ©. Although this 
coincidence is quite remarkable, it can not be treated as 
confirmation of the presented model, because it is not 
obvious that two dimensional MHD properly reproduces 
local structure of the three dimensional MHD turbulence 
in the place perpendicular to the mean magnetic field. It 
also should be noted that the results obtained in Ref. [34[ 
for the inertial range of two dimensional MHD turbulence 
are closer to the Kolmogorov spectrum. 

On the other hand, two dimensional MHD represents 
the system with two degrees of freedom (see, e.g., Eqs. 
(l)-(2) of Ref. [33j]) which keep several important proper- 
ties of MHD. Namely, linear modes have the same disper- 
sion and Alfven effect is also presented in two dimensional 
MHD. In this context theoretical as well as numerical 
study (including the case with strong background mag- 
netic field) of two dimensional MHD seems to have its 
own value, especially for the study of the turbulence lo- 
cality. 



IX. CONCLUSIONS 

Novel model of incompressible magnetohydrodynamic 
turbulence in the presence of a strong external magnetic 
field is proposed. It is suggested that in the presence of 
the strong external magnetic field incompressible MHD 
turbulence becomes nonlocal in the sense that the low 
frequency modes cause decorrelation of interacting high 
frequency modes from the inertial interval. Obtained 
nonlocal spectrum of the inertial range of incompress- 
ible MHD turbulence, given by Eqs. §8§ and (fTP)) . rep- 
resents the anisotropic analogue of Kraichnan's nonlocal 
spectrum of isotropic hydrodynamic turbulence. Based 
on the analysis performed in the framework of WCA it is 
shown that in contrast to the isotropic hydrodynamic tur- 
bulence, incompressible MHD turbulence could be both 
local and nonlocal and therefore an anisotropic analogues 
of both Kolmogorov and Kraichnan spectra are realizable 
in incompressible MHD turbulence. 
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APPENDIX A 

DIA equation that governs the dynamics of the Green 
function G h (k, t) for the inertial range of the isotropic 
hydrodynamic turbulence has the form (23l . [27j 

d f°° 

jG h {k, t-t') = - J d 3 k 1 d 3 k 2 5 k _ kl _ k2 6 1>2 x 
J dt"G h (k 1 ,t - t")Q h (k 2 ,t - t")G h (k, t" - t').(Al) 

The standard method of the analysis (2?| implies inte- 
gration of this equation with respect to r = t — t' in the 



interval (0, oo). Taking also into account Eqs. (|52|) - ([53)) . 
after the straightforward manipulations this leads to Eq. 
(|58| . Described method of derivation of Eq. (|58|l is not 
unique. Before performing integration with respect to r 
we can multiply both sides of Eq. (|A1[) by any function 
of t. For instance, if we multiply Eq. (IA1|) by exp (—at) 
and then repeat the same procedures, we obtain Eq. (|57p 
with uj replaced by ia. The reason why this equation is 
'wrong', whereas Eq. (|58f is 'right' is the same as in 
WCA analysis presented in Sec. [V] it is implied that r ac 
is the shortest timescale presented in the problem which 
totally determines the dynamics of the turbulence. 
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